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We re-examine the electrodynamic Casimir effect in a wedge defined by two perfect conductors 
making dihedral angle a — it /p. This system is analogous to the system defined by a cosmic string. 
We consider the wedge region as filled with an azimuthally symmetric material, with permittiv- 
ity/permeability ei,fii for distance from the axis r < a, and £2, 1^2 for r > a. The results are closely 
related to those for a circular-cylindrical geometry, but with non-integer azimuthal quantum number 
mp. Apart from a zero-mode divergence, which may be removed by choosing periodic boundary 
conditions on the wedge, and may be made finite if dispersion is included, we obtain finite results 
for the free energy corresponding to changes in a for the case when the speed of light is the same 
inside and outside the radius a, and for weak coupling, |ei —£2] <C 1, for purely dielectric media. We 
also consider the radiation produced by the sudden appearance of an infinite cosmic string, situated 
along the cusp line of the pre-existing wedge. 

PACS numbers: 42.50.Pq, 42.50.Lc, ll.lO.Gh, 98.80.Cq 



I. INTRODUCTION 

Quantum field theory in the wedge geometry contin- 
ues to attract interest, especially in connection with the 
Casimir effect. Usually it is assumed that the interior 
region of the wedge is a vacuum, and that the two plane 
surfaces 6* = and 9 = a {a denotes the opening angle) 
are perfectly conducting. The coordinate system is con- 
ventionally oriented such that the z axis coincides with 
the the singularity axis, i.e., the intersection line for the 
planes. For an introduction to the wedge model one may 
consult the book of Mostepanenko and Trunov [l| . 

The Casimir energy and stress in a wedge geometry was 
approached already in the 1970s Since that time, 

various embodiments of the wedge with perfectly con- 
ducting walls have been treated by Brevik and co-workers 
0, IE and others 0] ■ More recently a wedge intercut 
by a cylindrical shell was considered by Nesterenko and 
collaborators, first for a semicircular wedge Q, then for 
arbitrary dihedral angle 0]. Local Casimir stresses were 
examined by Saharian and co-workers [13, [HI [13 • Rosa 
and collaborators studied the interaction of an atom with 
a wedge [H, [l3| , the situation under which the closely re- 
lated Casimir-Polder force was investigated by Sukenik et 
al. some years ago 15]. That interaction was first worked 
out by Barton jl6i]. 

One reason for the interest in the wedge geometry is 
the similarity with the formalism encountered in Casimir 
theory of systems having circular symmetry. This applies 
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to the case of a perfectly conducting circular boundary 

[ia[ii,[ii,[2i, as well as to the case of a dielectric cir- 
cular boundary [U [H [H, [H, [H] . Another reason for 
studying the wedge is the analogy — at least in a formal 
sense — with the theory of a cosmic string (cf., for in- 
stance, Ref. [IHl, or Ref. 0]). Let us briefly elaborate 
on the last-mentioned point. The line element outside a 
cosmic string is, in standard notation. 



-dt^ + dr^ + (1 - AGMfr^de^ + dz"^, (1.1) 



where G is the gravitational constant and M the string 
mass per unit length. This is the geometry of locally flat 
space, with a deficit angle <& = SttGM being removed. 
Let us introduce the symbols (3 and p by 
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Now comparing the electromagnetic energy-momentum 
tensor outside the string [23| 
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with the electromagnetic energy-momentum in the wedge 
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we see that /3 corresponds to p. Hence the deficit angle <& 
corresponds to 27r — 2q;. We shall return to this analogy 
later. Note that the stress tensor diverges at r = 0, which 
makes the definition of a total Casimir energy in these 
configurations problematic. (Possible solutions to this 
problem were offered by Khusnutdinov and Bordag [28|.) 
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FIG. 1; The geometry considered in sections Ulland lllll There 
is a cylindrical perfectly conducting shell at radius a. In these 
sections the indices of refraction are equal, = eini = t2fJ,2- 

A particular variant of the wedge model occurs if we 
introduce a cylindrical boundary of radius a in the cavity. 
The situation is sketched in Fig. [1] This model has been 
studied in particular by Nesterenko et al. and by Saharian 
et al.; cf. Refs. B i, [O, d with a wealth of further 
references therein. (For example, the fermionic situation 
for the circular case was discussed by Bezzera de Mello 
et al. [1^.) The model can be looked upon as being in- 
termediate between that of a conventional wedge, and an 
optical fiber. And that brings us to the main theme of the 
present paper, namely to study the situation of Fig. [1] in 
the presence of a dielectric medium, both in the interior 
r < a as well as in the exterior, r > a. We designate the 
two regions by indices 1 and 2. Thus in the interior the 
refractive index is ni ~ ^ei/ii with ei and fii being the 
permittivity and the permeability, whereas in the exterior 
we have analogously n2 — ^ti[i2- We take all material 
quantities ei,/ii and e2,/i2 to be constant and nondis- 
persive. The special case when the circular boundary is 
perfectly conducting is included in the general situation 
when there is simply a dielectric/diamagnetic boundary. 
The plane surfaces d — ^ and 9 = a are taken to be 
perfectly conducting, as usual. 

We begin in the next section by considering the Fourier 
decomposition of the TE and TM modes when the cir- 
cular boundary is perfectly conducting. This is the sim- 



plest case. Then we move on to give an expression for 
the Casimir energy. The case of a dielectric / diamagnetic 
boundary is considered thereafter. The results for the 
wedge are in general divergent, not because of the di- 
vergence associated with the apex of the wedge, which 
does not contribute to the outward stress on the circu- 
lar arc, but because of the corners where the arc meets 
the sides of the wedge. This divergence may be isolated 
in the azimuthal zero modes, independent of the angu- 
lar coordinates. We propose isolation and removal of 
this divergence; alternatively, if the perfectly conducting 
boundaries aX 9 — Q,a are replaced by periodic bound- 
ary conditions, these divergences disappear. When ei- 
ther of these devices are employed, we obtain numerical 
results for the resulting finite Casimir energy, referring 
to the boundary between the two regions, r < a and 
r > a, both for weak and strong coupling. Finally, we 
exploit the analogy with a cosmic string to calculate, via 
the Bogoliubov transformation, the production of elec- 
tromagnetic energy associated with a "sudden" creation 
of the full wedge situation, as compared with the initial 
case of a single-medium filled wedge. 



II. ZERO-POINT ENERGY IN THE INTERIOR 
REGION— PERFECTLY CONDUCTING ARC 

As mentioned, we consider an isotropic and homoge- 
neous medium with permittivity ei and permeability /ii 
enclosed within a wedge region limited by the conducting 
plane surfaces 9 = and 9 = a {< 27r). In an xy plane, 
the cusp is situated at the origin. We use cylindrical co- 
ordinates (r, 6*, z). We employ Heaviside-Lorentz units, 
and put h and c equal to unity. 

Assume, to begin with, that the wedge is closed by 
a perfectly conducting singular arc at r = a. We write 
down the fundamental modes for stationary electromag- 
netic modes in the interior wedge, by invoking the ex- 
pansions given in Ref. (sdj : 
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(2.1b) 
(2.1c) 
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(2.1e) 
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Hz =2 ^ Jmp(Air)6'„Fo cosmpO. 

m=0 

I 



(2.1f) 



Here k is the axial wave number, and Ai is the transverse 
wave number given by 



A? 



2 2 



(2.2) 



The Jmp's are ordinary Bessel functions of order nip, 
which are finite at the origin for mp > 0, for p non- 
integral, while 



Fq ~ exp(i/cz — iujt) 



(2.3) 



is the m — version of the more general quantity Fm ~ 
exp(imp9 + ikz — iujt). The expressions (|2.1ap - (|2.1f1) 
satisfy the electromagnetic boundary conditions on the 
surfaces 9 = and 9 = a automatically, for arbitrary 
values of the coefficients am and bm- The i superscript 
on the coefficient refers to the interior region. The 
modes and the bm modes are independent of each other. 

Because of the closure of the region at r = a the prob- 
lem becomes an eigenvalue problem. Only discrete values 
of the transverse wave number Ai can occur. Let us dis- 
tinguish between the two kinds of modes: 

(i) TM polarization (the a5„-modes), which correspond 

to 



For m = only the TE mode is to be included. We 
now make use of the argument principle. Any meromor- 
phic function g{uj) of a complex variable w satisfies the 
equation 

where Wzoros are the zeros and ojpoics the poles of g{uj) 
lying inside the integration contour. The contour is cho- 
sen to be a large semi-circle in the right half-plane with 
radius R, closed by a straight line along the imaginary z 
axis from lu = iR to uj = —iR. A general advantage of 
this method is that the multiplicities of zeros as well as 
for poles are automatically taken care of. 

In the present case it is evident that g{u)) can be chosen 
as the product of Jmp and J^p. There are no poles in- 
volved, and the contribution of the large semi-circle goes 
to zero when R ^ oo. Thus we obtain 
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X —ln[Jrnp{Xia)J'mp{Xia)] . (2.10) 



Jmpi^ia) — 0. 



(2.4) 



We denote the roots by jmp,s, where s = 1,2,3.... For a 
given value of the axial wave number k the energy eigen- 
values are accordingly 



TM 

msk 



1 

nia 



Jmp^s 



k^a^, m>l, s>l. (2.5) 



(ii) TE polarization (the 6*„-modes), which correspond 



to zeroes of J,'„p, 
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nia V ^ 



m>0, s>l. (2.6) 



The interior zero-point energy per unit length is 
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(2.7) 

We here include the zero-point energy associated with 
the azimuthally symmetric TE mode, although there is 
no such TM mode. 

To simplify the formalism somewhat, we introduce the 
symbol 



In the second integral, k and lo are to be regarded as in- 
dependent variables in Ai = Ai(fc,Ci;). We now introduce 
the imaginary frequency C via lo — > i^, whereby 



Ai = y^nV _ fc2 ^ y-{nlC + k^) = (2.11) 

We thus encounter Bessel functions of imaginary argu- 
ments, Jmp{ix), with X = Kia, m > 0. Introducing the 
modified Bessel function I,^{x) via J^{ix) = i'^Iv{x) for 
arbitrary order v we get 



£. 



int 



c^CC;77ln [lmp{x)I'mp{x)\ 



l]_ f°° dk 

2 2^7_^2^„ ^ 

(2.12) 

Here we rewrite the derivative as d/d( = [n\a^C^/x)d/dx, 
and take into account the symmetry properties 
dk 2 dk, dC - 2 dC to get 



27r2 



dk 



edc 



^ 

X ^\n[lmp{x)Cp{x)] . (2.13) 

In the plane spanned by the axes niC, and k we may 
introduce polar coordinates X and Y 
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X = niC = Ki cos ( 
Y = k = Ki sin 9, 



(2.14a) 
(2.14b) 
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fulfilling the relation X"^ + = k\. The area element 
in the XY plane becomes KidnidO — riidC^dk. The in- 
tegration of the polar angle over the first quadrant then 
becomes simple, /J^^^ cos^ OdO — it/A, and we get 
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1 



dx 



dx x^ 



[-^nip('^)-^mp('^)] 
Imp{x) I'{x) 



(2.15) 



Going back to Eq. (|2.7p we can thus write the interior 
zero-point energy as 
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dx X 



(2.16) 



where the last term represents the TE m ~ mode. No 
regularization procedure has been applied at this stage. 



III. EXTERIOR REGION INCLUDED, 
ASSUMING PERFECTLY CONDUCTING 
CIRCULAR ARC 

We now include the exterior region r > a, still assum- 
ing the circular arc at r = a to be perfectly conducting. 



A choice has to be made for what kind of medium to 
fill the space r > a. One possible choice might be to 
assume a vacuum on the outside. Another natural choice 
would be to take the exterior medium to be identical 
to the interior one. We will in this section allow for a 
generalization of the last option, namely to assume that 
the exterior space is filled with a medium with arbitrary 
constants £2 and fi2, but with the restriction that their 
product is the same as in the interior: 



(3.1) 



We will refer to this situation as "diaphanous." This con- 
dition implying the constancy of light everywhere has un- 
der several occasions turned out to be convenient math- 
ematically, for instance in connection with the Casimir 
theory for dielectric balls [Ml, [H [H H HE HE HI, 
and in the Casimir theory for the relativistic piecewise 
uniform string [H, [s^, [io, ST], H^] (a review is given in 
Ref. Q). In the latter case, the velocity of light is to be 
replaced with the velocity of sound. The condition (|3.ip 
means in the present problem that A takes the same value 
on the outside as on the inside (assuming k to take the 
same values on the two sides). The principal advantage 
of this assumption, which is not easily satisfied in nature, 
is that in simple cases Casimir self-energies will turn out 
then to be finite. 

In the exterior region r > a we have the expansions, 
keeping the formalism at first quite general. 
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Fq cos mp9, 
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Fq cos mp9, 
Fq sin mpO, 



=2 Y,Hi]liX2rK,Fo cos mp9. 
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(3.2a) 
(3.2b) 
(3.2c) 
(3. 2d) 
(3.2e) 
(3.2f) 



As before, Fq is given by Eq. (|2.3|) . The presence of the 

Hankel function of the first kind, Hmp, ensures proper 
behavior (outgoing waves) at infinity. The e superscript 
refers to exterior modes. 

Let us now take into account the condition (|3.ip . im- 



plying Ai = A2 = A, and consider the boundary condi- 
tions. For the TM polarization (the a^-modes) we get 



i7«(Aa) = 0, m>l, 



(3.3) 
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whereas for the TE polarization (the 6^-niodes), 



0, m > 0. 



(3.4) 



The roots of these eigenvalue equations are complex — 
nevertheless, the argument principle may be applied as 
has been explained in detail in many places [ij, El, EH] . 
We can now calculate the exterior zero-point energy £°^^ 
in the same way as above. The modified Bessel function 
is introduced via Hl^\ix) = {2/T:)i-^''+^'^ K^{x). For 



the total zero-point energy/length E = 5'"* 
obtain 
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We now must face up to the fact that our result con- 
tains an unremovable divergence, associated with the 
nonzero 02 heat kernel coefficient found by Nesterenko 
et al. jl, Q . This occurs precisely because of the m = 
terms in Eq. (|3.5[) . If we were to write that expression as 
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terms are 



dx x^-' —\n\Ux)Kn{x)] 
dx 
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(3.8b) 



so the \/x terms cancel between the two modes (alter- 
natively those terms may be removed by contact terms, 
as we will see in the following), but the subleading 1/x'^ 
terms constitute an unremovable logarthmic divergence. 
Here, we have indicated an analytic regularization by tak- 
ing s to zero through positive values, which corresponds 
to the following divergent terms as s — > 0: 

I r°° d 1 

--/ dxx^-'—\nh{x)KQ{x) ^ —, (3.9a) 
z Jo dx as 

^ f°° d 3 

2], dxx^~^-\nI',{x)K',{x) ^ -. (3.9b) 

These precisely correspond to the two mode contribu- 
tions, adding up to l/2s, found by Nesterenko et al. Q. 

This zero-mode divergence is due to the sharp corners 
where the arc meets the wedge. We will proceed by set- 
ting this term aside, and computing the balance of the 
Casimir free energy. We note there is a closely related 
problem which Nesterenko et al. 0] dubbed a cone. That 
is, we identify the two wedge boundaries at = and 
a, and impose periodic boundary conditions there. This 
means that we may take the angular function in the mode 
sums to be e'™^^ , where m may be either positive or neg- 
ative, and where now p = 27: / a. Now all modes, includ- 
ing the zero modes (m — 0) contribute equally, and the 
summation on m becomes 



where the prime on the summation sign means that the 
m = terms are counted with half weight, we see in the 
following that the summation, £, may now be rendered 
finite (see Appendix A), but the residual correction, 5, 
is divergent. 

It is instructive to break up this residual zero mode 
contribution into its Dirichlet (TM) and Neumann (TE) 
parts. The former involves, asympotically for large x 



h{x)Ko{x) 
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2x 



while the latter requires 
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(3.7a) 



(3.7b) 



Then the two contributions to the residual zero-mode 



E 



(3.10) 



m— — 00 



m=0 



with the zero modes both having 1/2 weight in the latter 
form. (For the radial function in the interior we can 
only use in order that the solution be finite at the 
origin.) Thus we get precisely 2£ [Eq. (|3.6p ] without 
the residual zero mode term f , and we have eliminated 
the unremovable logarithmic divergence. This is because 
the sharp corners, where the arc meets the wedge, have 
been removed because there is no wedge boundary. So 
if the reader prefers, he or she may regard the rest of 
the discussion in this and the following section to refer to 
this situation, which will introduce an additional factor of 
two into the Casimir free energy, and with the restriction 
p > 1, where p = \ cor resp onds to the circular cylinder 
first considered in Ref . [l3| . 

So in any case disregarding in the following the residual 
zero- mode pieces £, we consider now the regularization 
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of the X]m=o' terms in Eq. (|3.6p . £, which, in order to 
be a Casimir energy, ought to be given in such a form 
that it reduces to zero in the limit when a ^ oo. This 
wiU chminate the divergence associated with the apex, 
which is not relevant to the force on the circular arc. It 
is easy to satisfy this requirement by observing that for 
large values of x, and for general v, we can approximate 

1 nr _ 

yznx V ^2; 

(3.11) 

implying that ^ I„ and ^ —Ki,. Accordingly, 
^ In i-I.KJlKl) ^ 2^ ln(/,A-,) ^ (3.12) 

to leading order in x. This term is to be subtracted off 
from the integrand in Eq. (|3.5p . The Casimir energy for 
the wedge becomes then 
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We may here perform a partial integration (the boundary 
terms at a; = and x = 00 do not contribute). 
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xdx 



X In [-Ax^I^p{x)I'^^{x)K^p{x)K'^^{x)\ . (3.14) 

It is helpful to introduce a quantity Xu{x) for arbitrary 
order 



K{x) = {h{x)K,{x))' , 



(3.15) 




FIG. 2: The wedge with a dielectric/diamagnetic boundary 
at r = a. In section ITVl we will allow ni 7^ 712. 



self-energy may be extracted from this formula. For now, 
we observe that this is a generalization of the standard 
formal result for a conducting circular cylinder, which is 
obtained from this result in the special case p = 1 [iTj . 
(The overall 1/n comes from an elementary scaling argu- 
ment (IBl-) When p = 1 the expression (|3.17p is one- half 
that for a conducting circular cylinder. Referring to the 
perfectly conducting wedge boundaries, we see that the 
Casimir energy for periodic boundary conditions, with 
period 2tt, is twice the Casimir energy for a perfectly 
conducting boundary condition imposed on the 9 inter- 
val of TT, a result obvious from the replacement of e""^ for 
m of either sign in the former case by smm9 or cos m9, 
m > 0, in the latter. This general observation, which 
is the theorem stated in Eq. (2.49) of Ref. [i^] (see also 
Ref. nil) will allow us to obtain numerical results rather 
immediately. For the periodic boundary condition situ- 
ation, which eliminates the zero-mode problem, p ~ 1 
is exactly the circular cylinder problem, and there is no 
additional factor of 1/2. 



IV. DIELECTRIC BOUNDARY AT r = a 



and to take into account the Wronskian W{I,j, K^} — 
— 1/x. From this we calculate the following useful rela- 
tionship 

- 4x^I,{x)Il{x)K,{x)Kl{x) = 1 - x^\l{x), (3.16) 

and so end up with the following convenient form for the 
Casimir energy: 
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00 
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: dx In [1 



x'Xlpix)] . (3.17) 



This is thus the boundary-induced contribution to the 
zero-point energy. If the boundary r = a were removed 
and either the interior or the exterior medium were cho- 
sen to fill the whole wedge region, we would get £ — 0. 
This is a property relying on the condition p.ip above. 
The temperature is assumed to be zero. 

Although the leading behavior of the Bessel functions 
has been subtracted in Eq. p.lTp , it is still not in general 
finite. We will see in the following section how a finite 



Assume now that the perfectly conducting arc at r = a 
is removed and replaced by a dielectric boundary, where- 
with the interior and exterior regions become coupled via 
electromagnetic boundary conditions at r = a. As before, 
we assume that the plane surfaces 9 — and 9 = a are 
perfectly conducting for all values of r. (Alternatively, 
we may impose periodic boundary conditions there.) 

We shall assume in the following that the media are 
arbitrary, with real and constant parameters ei, /ii in the 
interior and €2, fJ-2 in the exterior, without any restriction 
imposed on their product. This will, however, result in 
general in a divergent Casimir self-energy. 

Let A2 be the transverse wave number in the exterior 
region. 



2, ,2 



(4.1) 



with ~ The basic expansions are Eqs. (|2.1ap - 

(|2.1f|) in the interior and Eqs. p.2ap - (|3.2f|) in the exterior. 

As for the boundary conditions at r = o, only the tan- 
gential field components have to be taken into account. 
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From the continuity of and we get respectively 

J„p(u)aj„ = Hj^l{v)a^^ (4.2) 

and 

JmpiuW^ = H(^l{v)b'^, (4.3) 



where we have defined 



M = AiO, 



Xya. 



(4.4) 



Prom the component Eg we get 
kmp 



and from the component iJ^ 



kmp 



2 •^mp(u)6™ = —H^l'{v)a^.^ 



mp 



kmp 



(4.5) 



(4.6) 



The two last equations mean that a superposition of the TM and TE waves is in general necessary to satisfy the 
boundary conditions. The exception is the axially symmetric case m = 0. The condition for solution of the set of 
linear equations is that the system determinant vanishes. Observing the relation 



2 2 / 2 2n 2 2 

u — V ~ [n-^ — n2)U) a 
which follows from Eqs. (|4.ip and (|4.4p . we obtain after some manipulations the condition 



u Jmpiu) V H^liv) 



U Jmp{u) 



2 212 

m p k 



1 



(4.7) 



(4.8) 



r 



This is essentially the same transcendental eigenvalue 
equation as found for a step-index optical fiber (cf., for 
instance, Refs. [SO*] or [49*1 ). In transmission problems, 
one is usually interested in calculating the discrete values 
of the propagation constant /c, assuming that the waveg- 
uide is fed with some frequency lo. Here our intention is 
different, namely to calculate the discrete values of uj on 
the basis of an input value for the continuous axial wave 
vector k. As we noted in the previous section, this dis- 
persion relation generalizes that for a circular cylinder, 
the special case p = 1. 

It may be noted that the roots of Eq. (|4.8p are both 
real and complex. Application of the argument principle 
to such a problem is discussed in Ref . (50j ■ 



A. n\ — n2 

The TE and TM modes decouple in the special case 
when ni = y^eiA'T = n2 = ^/eiT^i- In this case, the 
dispersion relation (|4.8p reduces to AA = 0, where A and 
A are the two factors on the left-hand side of Eq. (|4.8p , 
and then using the Wronskian, we find after Euclidean 
rotation, u — > iC, 



AA 



iei+e2)l-ex^iImpK^pY^ 



4c2 



£162 



mp 



(4.9) 



where x = ko, c = l/n, and the reflection coefficient (for 
either polarization) 



£2 - £l 
£2 + ei 



(4.10) 



We conclude that the formula for the (zero-mode sub- 
tracted) Casimir energy is 



11 dk ^ , r'°° , d , , , , 



where 



gUx) = l-ex^>Lp. (4.12) 



where Amp is given by Eq. (j3.15p . Here, we have again 
subtracted off the terms that would be present if either 
medium filled the entire wedge. (The divergence struc- 
ture of the zero-mode term subtracted from Eq. (|4.1ip 
is analyzed in Appendix B.) Again cavalierly integrating 
by parts, we obtain, using the change of variables (I2.14p . 

^~-W^7 dxxHl-ex'Kap]- (4.13) 

m=0 •'^ 



As expected, this differs from the conducting case p.l7p 
by the appearance of in front of Amp- The conducting 
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case is obtained by setting ^ = 1. All of this is just as 
for the circular cylinder case, which is obtained from the 
p — I result by multiplying by a factor of 2. 

Let us now extract both the = 1 (perfect conducting) 
and the small ^ results for arbitrary p. A simple route 
is to follow the method given in in 18] or in Chap. 7 of 
Ref. j47j . The point is simply that the uniform asymp- 
totic expansion for the modified Bessel functions yields 
an asymptotic expansion for large p. Thus we can write 
(see the Appendix A for details) 



2n£ = 



e 



■ ln(27r/p) + f + 2 ^ £„ 



(4.14) 



where 



So = 



dx X 



ln(l-eVA2(x)) 



+ 



X 



f — 



47ra"^ 



4 (1 + 

oo 

dx X 



(4.15a) 







ln(l-eVA^p(x)) 



+ 



4 irn^p^ + x^Y" 



(4.15b) 



(Further details are given in the cited references.) Be- 
cause of the subtractions in the integrals, they are con- 
vergent. 

Let us first consider ^ as small, and keep only the terms 
of order Using the uniform asymptotic approximants, 
we find for large mp, 



1 



47ra2 V96toV 3840m V 
while numerical integration gives 



^0 = 



(-0.4908775). 



(4.16) 



(4.17) 



Thus 



-0.4908775+ - In 27r/p 



111 

288 p2 



77^4 



1 



M 



172800 p4 



"^^[fimp] - g{mp)] 



e 



re(p), 



(4.18) 



where we have added and subtracted the first two terms 
in the uniform asymptotic expansion. 



.9(^) 



1 



96i^2 38401/4' 
and / is the integral appearing in Em'- 



(4.19) 



0.0 



-0.2 



-0.4 




P 



FIG. 3: Casimir energy for weak coupling, <g 1, as a func- 
tion of p which is related to the dihedral angle a — n/p. (a) 
p > 1. (b) 0.5 < p < 1; the upper curve shows the exact 
energy, the lower the leading asymptotic approximation, ob- 
tained from Eq. (|4?T8)) by setting M = 0. 



In principle we are to take the Af — > oo limit in Eq. (I4.18P . 
In practice, we may keep only a few terms in the m sum. 
For example, keeping none of those corrections, that is 
setting M = 0, we get for p = 1, e(l) « -0.0010847. 
Keeping three terms is sufficient to find that e(l) is less 
than 10~^; indeed, the circular cylinder value is e(l) = 
[m m m H [23. This function e{p) is plotted in 
Fig. [3^, for p > 1, where it is sufficient to keep the 
leading asymptotic approximations; for p between 1/2 
and 1 (a between tt and 27r) we must retain at least 
one correction, M = 1, as shown in Fig. (No ob- 
servable change occurs with larger M.) Numerically, we 
see that the value for a cylinder with a conducting sep- 
tum {p = 1/2) is indistinguishable from e(0.5) = 1/4. 
Recall for periodic boundary conditions on the wedge 
(the "cone") p = 27r/a > 1, and an additional factor 
of two appears in the energy. 

Similarly, following the same references, we can obtain 
the strong coupling (perfect conductor) limit, ^ = 1. This 
time the formula for the energy is 



dx x^ 



4 (x2 + iy2)3 



(4.20) 



1 

87rna2 



e(p), 



(4.21) 
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where 



e{p) = 



1 7tT^ 1 TT^ 1 

-0.651752 + - In 2n/p H - 

4 '-^ 2880 352 32256 p4 

AI 

+ 2^ [/(TOp)-g(TOp)], (4.22) 

m— 1 

oo is understood. Now / is 



where again the Uniit M 
given by 



dx i 



1 



4 (x2 + jy2)3 



, (4.23) 



and now the asymptotic terms are 

7 5 



960i/2 



3584:/4 



(4.24) 



Keeping no correction terms is already very good at 
p = 1, where with Af = e(l)/(47r) « -0.013633, oiih^ 
shghtly different from the exact answer of —0.01356 [17l |. 
Keeping just AI — I gives exact coincidence to the indi- 
cated accuracy. This function e{p) is plotted in Fig. 0^ 
for p > where the asymptotic approximation is suffi- 
cient, while two correction terms are included in the re- 
gion 0.5 < p < 1, as shown in Fig. IHd. It is curious that 
the energy vanishes now not at p — 1, but at p — 0.583. 
(Again, recall only p > 1 is relevant for periodic bound- 
ary conditions on the wedge.) 



B. ni / n2, /ii = = 1 

Finally, we can follow Ref. [2l[ to obtain the weak- 
coupling Casimir sclf-cncrgy for a purely dielectric wedge, 
where /ii = /i2 = 1. We can only examine the coefficient 
of (ei — £2)^ because the result is divergent in higher 
orders. It is hardly necessary to give details, since all 
that is necessary is to replace m by mp in the analysis 
given in that reference. The energy per area in the wedge 
is 



2 °o r°° 



(ei - e2) 
327rna2 ^ 

m=0 



(4.25) 



where the exact form of gm is elaborate, but has the 
asymptotic form 



11 

9miy)- T;z:j:^Y.7ZZAkfkiz^' mp^oo, (4.26) 



2m? p^ ^-^ {mpY 
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FIG. 4: Casimir energy for strong coupling, = 1, as a 
function of p, related to the dihedral angle by q = n/p. (a): 
p > 1; in this region M = in Eq. (I4.22|) is sufficient, (b): 
0.5 < p < 1; in this region M = 2 in Eq. (|4.22p is sufficient, 
and the comparison with the M = Q result (lower curve) is 
made. 



where y — mpz, and are rational functions of z, given 
in Ref. [21], about which all we need to know here is 



p2 hm > 



,2-s 



m— 1 



m " / dz z ^/i(z) 




i C(3) 
167r2 



dz z f2{z) 



(4.27a) 
(4.27b) 



1™ 



m=0 



— V — / dzz^fsiz) 
P ^1 Jo 



dzz^-'h{z) = ^ln27r/p, 

(4.27c) 

dzz^f4,{z) = 0, (4.27d) 

197r2 1 , 
7680 p'' 



dzz^fQ{x) = 0, 



(4.27f) 



1^1 



dz z^friz) 



2097r'* 1 
' 5806080 p4' 

(4.27g) 



Here a contact term, which cannot contribute to any 
observable force, has been removed from Eq. (j4.27bp . 
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1.0 



FIG. 5: Casimir energy for weak coupling, |ei — £2] <C 1, as a 
function of p, which is related to the dihedral angle a = ir/p, 
for a purely dielectric wedge, (a): p > 1; only the leading 
asymptotic terms have been included, (b): 0.5 < p < 1, the 
effect of the remainder is significant; the upper curve shows 
only the leading asymptotic terms, while the lower curve in- 
cludes the remainder function r. 



Again, for the precise definition of Eq. (|4.27cp see Ap- 
pendix A. Then the Casimir energy per unit length of 
the dilute dielectric wedge is 



£ — —w{p), 



where 

w(jp) 

where 



C(3) 5 
167r2 32 



ln(27r/p) 



r{mp) - 



1971^ 

7680J52 

0.000012 

9 ' 



(4.28) 



- 0.301590 



(4.29) 



r{v) = 2 / dyy 



k=l 

(4.30) 

and we have used the next term in the asymptotic series 
to estimate the contribution for to > 5. This, numeri- 
cally, yields the correct value of zero for p = \. The values 
for the Casimir energy for larger values of p are shown 
in Fig. [5^, and for smaller values of p in Fig. ^p. For 
the septum case, the numerical value of w(0.5) = 0.1666, 
which seems likely to represent exactly 1/6. (Periodic 
boundary conditions restrict p > 1.) 



V. ENERGY PRODUCTION IN THE SUDDEN 
FORMATION OF A COSMIC STRING 

As already mentioned, the electromagnetic theory of 
the wedge is related to the theory of cosmic strings. In 
general, cosmic strings are believed to be possible in- 
gredients in the very early Universe; they are related to 
phase transitions. One particular aspect of this study is 
to estimate the energy production in the form of mass- 
less particles when a string is formed "suddenly" at some 
instant t = to, where to is a characteristic time usually 
taken to be of order 10"^*^ s as is typical for grand unified 
theories (GUTs). One can calculate the number of parti- 
cles associated with the formation of the string in terms 
of Bogoliubov coefficients relating the initial Minkowski 
metric to the (static) metric after the string has been 
formed. This approach was pioneered by Parker [sij for 
massless scalar fields, assuming the string radius to be 
zero. Analogous calculations were made in Ref. (s^ in 
the electromagnetic case, still assuming the string radius 
to be zero, and in Ref . 1531 for the case of a finite string 
radius. (See also Ref. [28j.) 

In the following we shall investigate the following 
model: Let a cosmic string of vanishing radius and large 
length L be formed suddenly along the cusp line of the 
wedge, i.e., along the z axis. The time of formation is 
t = to- In the interior region we assume that there is 
one single isotropic medium present, with refractive in- 
dex n = ^JTp,. The interior region is closed by a perfectly 
conducting arc with radius r = a. In the present context, 
a plays the role of a "large" boundary; in connection with 
strings, a is usually taken to be of the same order as L 
[Slj . The finiteness of a will moreover make it possible 
to normalize the fundamental modes in a conventional 
way. As for the description of the electromagnetic fields, 
we have to distinguish between the Minkowski metric 
present for t < to and the string metric for t > to- 



A. The case t < to 

We shall consider the TM mode only, for which the 
central field component is Ez- For definiteness we repro- 
duce the fields in the to, k mode here, in a convenient 
notation (replacing the previous combination 2ia^ with 
the symbol TV): 



= NJmp{Xr)Fo siiimp9, 



Er = —NJ!^p{Xr)Fosinmp9, 
A 



(5.1) 



(5.2) 



Ee = ——NJ^p{Xr)Fo cos mp9, (5.3) 



Hz = 0, 



(5.4) 
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impeuj 



N J,np{Xr) Fq cos mp9, (5.5) 



is the PoyntinK vector. As in any linear wave theory we 
can set [54| 



He = —NJ'^p{Xr)Fo simnpd. 



(5.6) 



As before, A = ^Ju^llP- — fc^, Fq — exp{ikz — iujt), 
and m > 1. The boundary condition on the arc is 
Jmp{Xa) = 0, giving the solutions Xms, s — 1,2,3,... 
for the transverse wave number A. 

It is now convenient as an intermediate step to make 
use of the formalism of scalar field theory. Define the 
scalar field mode 'ipmsk, satisfying Dirichlet boundary 
conditions on all surfaces, as 



Ipmsk = NJmp{Xmsr)Fo smmpO. 



(5.7) 



It is seen to have the same form as the m, s, k mode of 
the field component Ez. For reasons to become clear 
later, we choose the magnitude |iV| of the normalization 
constant N to be 



\N\ 



1 



Xr 



n \ aeujmsk a\Jmp+i{Xmsa)\' 



(5.8) 



with UJrnsk = (l/n)^/A^^ + fc2. 

We define the Klein- Gordon product as 



{'4^msk: ^m' s'k' 



A„ 



V'msfe 9o i/'m's'fe' rdrdOdz, 

(5.9) 

and then get by direct calculation 

{ipmsk ^i'm's'k') = 2'K6{k - k')Sjnm'Sss' ■ (5.10) 

Consider now the electromagnetic energy W in the wedge 
region. We may calculate it by integrating the energy 
density w over the volume: 

W ^ J wdV =^J [e\E\^ + fi\H\^]rdrd0dz, (5.11) 

using the general recursion equation (J^ = Ji,[x)): 



-^. = + (5.12) 



as well as the integral formula 

{Jmp-x{Xmsr) + Jmp^x{^msr)\ rdr 

= a J„^p^i{XmsO') , 



(5.13) 



which holds when Jmpi^msO-) ~ 0. It is however simpler 
to go via the axial energy flux P, given as 



P = / SzdA, 



(5.14) 



where dA — rdrdO is the cross-sectional area element, 
and where 



Sz = -^{ErHg — EeH*) 



(5.15) 



(5.16) 



where Cg is the axial group velocity. From Eqs. (j5.14p 
and (|5.15p we then get 



P = ^T^ 1^1 Jmp+li>^msa). 

OA™ a 



(5.17) 



In geometric units P has the dimension cm ^. As Cg = 
duj/dk = k/{n'^uj) we get for the energy per unit length 



W 

T 



2 2 9 

aen a uj^ 
8A^„ 



WJmp+li^msa). 



(5.18) 



We see that W/L is expressible in terms of |i?2p inte- 
grated over the cross section: 



W 

T 



771 Sk 



2X1 



\Ez\'dA. 



(5.19) 



This relation will turn out to be useful in the following. 

Quantum theory. We assume henceforth the real repre- 
sentation for the fields. The component Ez{r, t) = Ez{x), 
considered quantum-mechanically as a Hermitian opera- 
tor, is expanded as 



Ezix) 



J '^J2Y1 [(^^sk1pmsk{x) + a|„3feVmsfc(2;) 



m—1 s 



(5.20) 

where Umsk and aj^^j, are annihilation and creation op- 
erators satisfying the commutation relations 

^ -27rS{k-k')S„,^,S,s'- (5.21) 



We now go back to the relation (|5.19p . and require that 
the total energy W associated with the to, k, s mode is 
equal to the occupation number {o^^sk^msk) times the 
photon energy ujmsk- 



msk 



\2 



{EDrdrdOdz = (aLfe« 



msk 



k- (5.22) 



Here we insert the expansion (|5.20p . Because of the or- 
thogonality of (|5.21[) , the various modes decouple so that 
the total energy is a sum over the mode energies. For 
the mode ipmsk: written in the form (j5.7p . we get from 
the condition above the expression for the normalization 
constant \N\ already given in Eq. (jS.Sp . If n = 1 and 
a — 2tt, we recover the expression given in Ref. [s^ . 



B. The case t > to. The Bogoliubov transformation 

After the sudden creation of the cosmic string along the 
cusp line (the z axis) at the instant t = to, we assume 
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that the string metric is static. All transient phenom- 
ena are intended to be taken care of via the use of the 
quantum mechanical sudden transformation below. We 
first have to establish the field expressions in the presence 
of the string metric. The central gravitational quantity 
appearing in the formalism will be 



formulas. The normalization constant \N^\ becomes 



n V aeuji^sk a\J^p+i{Ai,sa)\ 
and the Klein-Gordon product, defined as 



(5.27) 



/3 = (1 - 4GA/)-\ 



(5.23) 



already introduced above in Eq. (jl.2b[) . In a string con- 
text, P is believed to be very close to unity. Writing the leads to 
field component as Ez{r) expiikz — icut) sinmpO we 
obtain the following equation for the quantity Ez{r): 



ipusk do V'^'s'fc' rdrdOdz, 

(5.28) 



1 d 



2 fS'^m^p'^ 



E,{r)=0, (5.24) 



^ dr^ r dr 

with = n^uj'^ — as before. Introducing the symbol 
v as 

V = (3m, 

we can write the fundamental v, s,k mode as 
i^iysk = Ni,J^p{X,,sr)Fo sinmp0, 



(ipiysk^ipi^'s'k') ^ 2TrS{k - k')S^^>dss'- (5.29) 
The quantum mechanical expansion for Ez becomes 

^J2Y1 [a-i'ski^i'skix) + aL/cCsfela;) , 

m=l s 

(5.30) 

(5.25) with associated commutation relations 

Ksk^<^l's'k'] = 27r(5(fc - k')S,,,6,,, . (5.31) 



(5.26) We have to specify the continuity conditions for the fields 
at the transition time tQ. The component Ez will be 
with Fq = exjp{ikz-iu!i,skt)- The boundary condition on required to be continuous, 
r = a is Jixp(Aa) = 0, giving solutions A^.^, s = 1,2,3... 
for the transverse wave number. 

The formalism now becomes quite similar to that given 
before in the non-gravitational case. We list the main 

I 



Ezix) ^Ezix) 



(5.32) 



as well as the Klein-Gordon product, 

/ Ezdo El rdrdOdz 

from which we get 



I3\ 



[Ez do El 



rdrdOdz, 



[doEz{x)],- = 



[doEz{x)\ 



(5.33) 



(5.34) 



The Bogoliubov transformation. We have now two kinds of basic modes, namely ipmsk for t < to, and ip^sk for t > to- 
There are correspondingly two vacuum states, satisfying the relations amsk\0)msk = and a,^sk\0)i/sk = 0. As in 
Refs. [52I, [5^ we may expand the modes in terms of each other: 



dk' 

^uskix) = / 7^ [jiiysk\m's'k')'ipjn>s'k'{x) + S{iysk\m's'k')i;l^,^,,,,{x)] , 

J-00 , , 

m's' 

where 7 and 5 are the Bogoliuobov coefficients [ssj . The corresponding expansions for the operators are 



(5.35) 



dk' 
2 



fc' r 

— ^ ^'y{vsk\m's'k')am's'k' + ^* {i^sk\m' s' k')al^, ^,f., 



(5.36) 



It means that the average number of particles produced in the m, s, k mode per unit k space interval becomes 



dk 



rik' 



(5.37) 
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From Eq. (|5.35p we obtain, when making use of the normahzation of the scalar product corresponding to string space, 



i'usk do 1pm' s'k' rdrdOdz. 



(5.38) 



Here we insert the expressions (|5.26l) and (|5.7p for ip^sk and ipm's'k', and for simphcity we put io = 0. Defining the 
quantity I^s' as 



Is 



Jo Jvp[Ksr)Jmp{>^ras'r)rdr 



we then obtain after some calculation 
6 {i' skim' s'k') = 



1 A,„, 



■ 2n6{k + k')S„ 



^vsk 
^msk 



^rask 



^ysk 



(5.39) 



(5.40) 



As the value of f3 is very close to unity, we put P — 1 
everywhere except in the difference between the square 



roots. With J^, 



Jmp and \i/s' 



Xms' , the numera- 



tor in (|5.39p reduces to (a^/2) J^p^j(Amsa)(5ss/, so that 
approximately 



(5.41) 



Moreover, by applying the integral operator J dk' jln on 
[27r(fc + fc')]^ we obtain effectively the length L of the 
string. For the electromagnetic energy produced in the 
mode m, s, /c, per unit wave number interval, we then get 

d^msk ^msk dNmsk 1 / 

dk L dk A \uJmsk i^i^sk 

(5.42) 

There are two properties of this expression worth notic- 
ing: 

(1) It is independent of the opening angle a. The phys- 
ical reason for this appears to be related to the fact that 
our region of quantization is the interior wedge region 
only. All the produced energy is taken to be channeled 
into the wedge region (we are thus not cutting out a 
fraction a/2TT of the total produced energy). This con- 
trasts the behavior in the cylindrically symmetric case, 
where the produced energy is azimuthally symmetric in 
the whole region < 9 < 2tt [sij. 

(2) The produced energy, when expressed in terms of 
frequencies, does not contain the refractive index n ex- 
plicitly. Equation (|5.42p is formally the same as Eq. (52) 
in Ref. [s^ . 

We may process the expression further by making use 
of the asymptotic formula for the roots of the Bessel func- 
tion. 



. , 1, TT 

X,nsa = STT + (to - -)-. 



(5.43) 



Here it is of physical interest to consider the region 
around zero axial wave number, k ~ 0. Then Lu^sk 



^lyso = Ks/n, ujmsk —> t^msO = Kns/n, leading to 

TO 



^msQ 



2s + TO 



(5.44) 



where we have expanded in the small quantity (/3 — 1) to 
second order. Then, 



dW, 



msk 



dk 



" 1 - (5.45) 
ona 2s + m — -k 



We thus see that finally the factor n turns up in the 
denominator; this is a characteristic property of Casimir 
energy expressions for dielectrics [46]. 
The simplest possibility m = s = 1 yields 



dWi 



Ik 



dk 



_ 47r ,.2 



fe«o 



.{P-XY^ (GUy. (5.46) 
zUna 5na 



The total energy W produced per unit length follows by 
multiplying (|5.46p with the wave number width A/c ~ 
1/L ^ 1/a around fc = 0. We may take a to be of 
the same order as the horizon size ^ t being the time 
just after the Big Bang. We thus get, when leaving n 
unspecified. 



(5.47) 



This is a characteristic property of cosmic string theory. 



VI. CONCLUSIONS 

We have computed the Casimir free energy for a wedge- 
shaped region bounded by perfectly conducting planes 
meeting in an angle. The wedge region is filled with 
an azimuthally symmetric material which is discontinous 
at a radius a from the intersection axis. In general 
the wedge geometry is plagued with divergence prob- 
lems. Familiar is the divergence associated with the apex, 
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which is not relevant to the force on the circular bound- 
ary. But there are also divergences associated with the 
corners where the circular arc meets the wedge bound- 
ary. These divergences are manifested only in the m = 
modes, which possess no dependence on the angular co- 
ordinate, and have here been isolated and disregarded 
in the calculational part of this paper. They will not 
be present if the perfectly conducting boundary condi- 
tions on the wedge are replaced by periodic boundary 
conditions, which restricts the parameter p to be greater 
than unity. Then, if the speed of light is the same both 
inside and outside the radius a, the energy correspond- 
ing to changes in a is finite. If the speed of light differs 
for r < a and r > a, the Casimir energy is finite only 
through second order in the discontinuity of the speed of 
light. These results are seen to be straightforward gen- 
eralizations of results holding for dielectric/diamagnetic 
circular cylinders, which are recovered if p = 1. We also 
consider, in the "sudden" approximation, the electromag- 
netic radiation produced by the appearance of a cosmic 
string in this geometry. 
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where we have let in the m terms x = mpz. Now the last 
factor is, as s ^ 0, 

l-|-2C(s)p"'' ^ -s(ln27r-lnp), (A4) 

while the integral diverges as s 0: 



dx 



1 



(1+^^)^ 



(A5) 



The result (|4.14p follows immediately. An identical argu- 
ment leads to Eq. (|4.27cp . This argument demonstrates 
the importance for achieving a finite result of including 
both TE and TM zero modes with half weight. 



APPENDIX B: CONVERGENCE CONDITION 
FOR ADDITIONAL ENERGY TERM ASSUMING 
HIGH FREQUENCY TRANSPARENCY 



The energy expression ()4.13p . 



1 ^ 

m=0 



has an additional term consisting of one half the m — Q 
term for the TE mode minus one half times that of 
the TM mode. These modes are determined in the di- 
aphanous case by A A — 0, where A and A are the two 
factors in Eq. (|4.8p . which for the zero modes are propor- 
tional to 



APPENDIX A: ANALYTIC REGULARIZATION 
OF LOGARITHMICALLY DIVERGENT TERM 



AoAoocfl^-l^^M-s.^V (B2) 



The only subtlety in the numerical calculations in 
Sec. IIVI is how the superficially logarithmically divergent 
terms are regulated. Starting from Eq. (|4.13p we have 



m=0 



where 



dx 



47ra2 7o 4(1 + 2:2)3' 



dx- 



(Al) 



(A2a) 



(A2b) 



IQ A{m?p'^ + a;2)3 

Here, it will be observed that the integrals over x are 
logarithmically divergent. We will regulate them analyt- 
ically by replacing in the numerator of both x^ — > x^~'^ , 
where we will at the end take s to zero through positive 
values. Thus we have 



dx- 



167ra2 Jo (l + a;2)3 
(l + 2f;(mp)-M ,(A3) 



Then, using the Wronskian, we see that the residual zero- 
mode term is 



1 



' 2 d , 1 + ^x\q{x) 
dxx -j-ln- . (B3) 

dx 1 — t,XAo{x) 



(This just says that the reflection coefficients for the two 
modes are ^tm = £, and ^te = —(,■) 

As in the perfectly conducting case, this is divergent, 
if ^ is constant, because 

{Io{x)Ko{x)y ^ x^^, (B4) 



which means that the integral in Eq. (jB3p is linearly di- 
vergent. However, if ^ is frequency dependent, so that 

^-r", C-oo, (B5) 
it is apparent that the integral becomes finite if /3 > 1. 
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